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THE SPANI7ISE LIFT DISTRIBUTION* 
By 0. Schrenk 

PRELIMINARY REMARKS 

In this paper a simple approximation method is pre- 
sented for rapidly computing the lift distributions of 
arbitrary airfoils. The numerical results are compared 
with those obtained by an exact method and for many pur- 
poses show a satisfactory degree of accuracy. The latter, 
for all practically occurring cases, can be estimatod at 
the start of the computation work with the aid of the com- 
parison examples given. 

The method described below enables the approximate 
determination of the lift distributions in a few minutes 
with an accuracy sufficient for many purposes. It is also 
characterized by a certain simplicity which is useful in 
the clarification of many questions and is in accord with 
the engineer's point of view. Finally, the method is ap- 
plicable to cases for which all other methods are entirely 
unsuitable (for example, wings with end plates). 

Similar methods, as the author subsequently found, 
have already appeared elsewhere. The surprising accuracy 
of such simple methods is, however, generally unrecognized, 
so that a presentation of comparison computations which 
provide a measure of the degree of accuracy obtainable, 
should lead to an extended application of the method. 

The author wishes to express his appreciation to his 
co-workor at Cottingen, Mr. N . Hiorth, for carrying out 
the laborious computations required for the comparison. 



*»Ein einfaches Naherungsverf ahren aur Ermittlung von Auf- 
triebsverteilungen langs der Tragf luigelspannweite . 
Luftwissen, Bd. 7, Nr. 4, April 1940, pp. 118-120. 
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FUNDAMENTAL IDEA 07 THE METHOD 



The plausible assumption is made that the real lift 
distribution lies "between an ideal distribution independ- 
ent of the wing shape and a distribution determined in a 
simple manner by the wing shape. The ideal distribution 
is that with minimum induced drag and constant induced 
downwash velocity - that is f for the usual monoplane, the 
elliptic distribution; while the distribution dependent 
on the the shape is proportional to at at each position 
of the wing* 

COMPUTATION PROCEDURE 



In the case of the untwisted wing the angle of attack 
is not absolutely necessary. There is drawn instead (for 
the monoplane) the semiellipse of equal area with the 
chord-distribution curve, and the lift distribution is ob- 
tained by forming the arithmetical mean between the two 
curves. 

In the case of twisted wings (and similarly for wings 
with aileron or flap deflection) there is first determined 
the zero lift direction of the entire wing. A sufficient 
approximation for this is the direction of the mean aero- 
dynamic twist 



j «(x) t(x) 

0 a — 

t 

Where the bars denote mean, values, 6(x) and "3 are 
the twist angles between the reference direction of the 
entire airplane and the local and mean zero lift direc- 
tions, respectively. For all further computations, the 
angles of attack and twist are reckoned from the zero lift 
direction given by 5. 

With angles reckoned as indicated above, the lift is 
decomposed into a component without twist and a twist com- 
ponent without lift. The second component is determined 
on the basis of a mean value which has the zero line in- 
stead of the ellipse as the ideal distribution, and for 
which the twist angle must always be considered. 

The general case with twist is most conveniently com- 
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puted with the aid of the following formula, which requires 
no explanation: 



dA _ 1 dc 
dx = 2 q da 



+ l q . *£» 8(x) t(x) 

In the ahove relation ~ a and ~ a can generally 

da da^, 

"be taken as constant along the span. A somewhat more accu- 
rate formula is given below. 

dc_ dc a 

The trial computation of and -g^- may he some- 

what refined "by the introduction of a correction factor 
(1 + k) setting : 

dc a 

dc a da o, 



da dc„ p 



x + - (i + K ) 



da^ ttId 



The constant K for various taper ratios was determined 
by trial in such a Tray that the lift determined with our 
approximation method is given as correctly as possible 
(fie;. !)• In the case of the rectangular wing the value 
thus determined agrees with the theoretically determined 
value of Glauort, For other taper ratios the agreement 
has not he en checked and in this connection is not required 
A further refinement in the value of K, nevertheless, 
seoms unnecessary • 

The computation procedure thus consists of tho follow- 
ing steps: 

a) Computation of 8 hy forming the mean. 

t>) Trial computation of dc a / da and dCg/do^. 

c) Computation of the lift distribution hy the 
formula for dA/dx. 
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ACCURACY OF THE METHOD 



The comparisons given in figures 2 to 13, "between the 
accurate values computed "by the method of Multhopp and the 
results of the approximation method, show in general a sat- 
isfactory, and to some extent even surprising, agreement. 
The error arising through the assumption: of a constant 
value of dc a /da along the span in unfavorable cases, can 
"be eliminated "by computing a mean value 



dc„ (x) 



dc 



a 



t(x) 



doQ, 



and then computing the lift "by the only slightly altered 
formula 



dA 1 

dS = 2 q ^ 



dc a<*> 

da 



da n J \| 



The sha,rp difference between lift without twi3t and twist 
without lift, to be sure, no longer arises. This correc- 
tion has not "been applied in the example here given* 

Further refinement through additional computations, 
at the expense of briefness and simplicity, does not ap- 
pear to he of advantage. The deviations, moreover, occur 
at such positions whero the exact theoretical solution 
does not agree with actuality.- for example, for unrounded 
wing tips and at transition positions of ailerons and 
flaps. 

The deviations at the transition positions of ailer- 
ons and flaps can ho readily balanced by hand \rith the aid 
of examples ?iven in the figures. 

. POSSIBILITIES OF APPLICATION 



The method, as may be seen, is suitable for all prob- 
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lems where too great accuracy is not required; that is, in 
general, for investigations with regard to the maximum 
lift coefficient c a , stalling, and static equilibrium 



' a max' 



problems. By the decomposition into an ideal, plan form, 
and twist distribution, simple and time-saving relations 
may be set up for frequently repented computations of bend- 
ing moments, transverse forces, and torsional moments along 
the span. For the computation of the induced drag and for 
the downwash computation, the method is not directly appli- 
cable. 

The method is suitable for the determination of lift 
distributions also in cases for which the usual methods 
fail completely. Thus, it is applicable to monoplane 
wings with end plates. Ideal distributions that take the 
place of the ellipse can be computed on the basis of an 
investigation of the Aerodynamic Experimental Institute 
(reference l) - the distributions there given being for 
smallest induced drag and constant downwash. The relation 
for the lift distribution now becomes: 

dA 1 dc a » . - fJx) 



a t(x) + t 
dx 2 da b 

2 



1~ / f (x) d(x) 

\b 

2 

+ | , ,(,) t(x) 

where f(x) is the function denoted by Mangier as the 
"ideal function" for the given case with end plates, 

She method should likewise find application to "bi- 
planes and other arrangements. 



SUMMARY 



The approximation method described makes possible 
lift-distribution computations? in a few minutes. Compar- 
ison with an exact method shows satisfactory agreement, 
The method is of greater applicability than the exact 
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method and includes also the important case of the wing 
with end plates. 



Translation hy S. Roiss, 
National Advisory Committee 
for Aeronautics. 



REFERENCE 



1. Mangier, W.: The Lift Distribution of Wings with End 
Platos. T.M. No. 856, NACA 1938. 



NACA Technical Memorandum No. 948 



Fig*. 1,2,3, 4, 5,6 
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Figure 1*- The values of K 

that occur in 
the formulas for dA/dx. For 

non-trcpazoid shaped wings 
the nearest value of the 
taper ratio is used for the 
determination of K. 
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Figures 2 to 5.- Co3??.rison curves with b 2 /F=6.67 for various 

t-.poT ratios. The comparison fo» ta.ti=0.5 is 
given on Fig. 6. TSo tA'ist, i 
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Figure 6.- Comparison curves for t a tti=0.5 for various" values 
of F/b 2 . No twist. 
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Pigs. 7,8, 9, 10 
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Figure 7.- Comparison curves 

for a wing without 
twist with cut-out in center. 
All three curves enclose 
equal area. 

Key for Pig. 7. 

Chord distribution 

Approximation 

Exact by Multhopp's 

method. 

Figure 8.- Comparison curves 
for a rectangular 
wing with semicircular 
rounding, without twist. The 
deviation of the approximation 
method at the wing tip for the 
non-rounded rectangular wing 
here almost completely 
disappears. 
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Figure 9.- Comparison curves 

for a wing with 
t a :ti=0.5 and b 2 /F=6.67 for 
various values of Ca with a 
linear twist which is 0° at 
the center and 3° downwards 
at the wing tip. 




Figure 10.- Larger scale 

representation 
of the curves for c a =0 frem 
Fig. 9. 
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Figs. 11,12,13 



Figure 11.- Results for c a =0 

as in Fig, 10 with 
a twist which similarly amounts 
to 3° at the wing tip but 
increases parabalically. 
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Figure 12,- Comparison curves for a rectangular wing with b /F=5 for 

an aileron deflection 0=0.25*14° with and without lift. 
The comers of the approximation curve can practically be well rounded 
off by hand following the example given and thuB considerably better 
agreement is obtained. The rolling moments of the non-rounded off 
approximation curve very well agree with the exactly computed values. 
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Figure 13.- Comparison 

curves for a 
tapered wing with t a tti= 
0.5 and b 2 /F=€.67 with 
60° deflected split flap 
along center half of 
span. 
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